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Abstract 


Processes which occur within the region between ^Rq and 25 Rq, 
which is called the solar envelope, probably have an important effect on 
the solar wind as seen at 1 AU. In the envelope the wind speed becomes 
supersonic and super- Alfvenic , the magnetic energy density is larger than 
the flow energy density, and the magnetic energy density is much larger 
than the thermal energy density. Large azimuthal gradients in the bulk 
speed are expected in the envelope, but the stream interactions near the 
outer edge of the envelope are probably relatively small. 

Cosmic ray observations suggest the presence of hydromagnetic waves 
in the envelope. The collisionless damping of such waves could heat protons 
out to «25 Rq and thereby cause an increase in V and consistent with the 
observed T^-V relation. A mechanism which couples protons and electrons 
would also heat and accelerate the wind. Alfven waves can accelerate 
the wind in the envelope without necessarily causing heating of protons; 
the Lorentz force might have a similar effect. 

Magnetic stresses cause an increase in which reaches a maximum in 
the envelope. Viscous stresses would also increase V , but their 
importance is controversial. 

Observations at 1 AU are not sufficient to determine the relative 
importance of the various processes that might occur in the envelope. 

Space probes which descend close to the envelope should significantly 
increase out understanding of such processes. 
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I. Introduction 

Several types of evidence suggest that a variety of physical processes 
which are important in the formation of the solar wind occur in the region 
between ^21^ and £=251^. This region, called the solar envelope, bridges 
the gap between the interplanetary medium and the lower corona (Figure 1) . 
Unlike the interplanetary medium, which has been studied directly by space 
probes to .7 AU, and unlike the lower corona, which has been extensively 
studied by various indirect methods, this intermediate region has remained 
largely inaccessible to observations. Yet it undoubtedly holds the answers 
to many of the unsolved problems in solar wind theory and solar-terrestrial 
relations, and it is probably the site of novel and interesting physical 
processes. 

The observations near 1 AU and at the sun provide valuable clues to 
the nature of the processes occurring near the sun, but they are not 
sufficiently restrictive to unambiguously identify the most important 
processes governing acceleration, heating, energy transport, and angular 
momentum transport in the envelope. In fact, there are competing theories 
for each of these processes. Without a satisfactory physical understanding 
of the processes one cannot derive an accurate model for the fluid parameters 
in the envelope. However, one can obtain an estimate of these 
parameters using models designed to explain the "quiet time" solar wind 
at 1 AU and from this one can obtain estimates of the parameters that 
govern the general physical processes in the envelope. This and the question 
of azimuthal and temporal variations in the envelope are the subjects of 
Section II. 

Section III contains a general discussion of the physical processes 
which might occur in the envelope. It begins with a discussion of cosmic 
ray propagation in the envelope, since the existence of an envelope was 



suggested by cosmic ray observations. These observations indicated a 
region of hydromagnetic wave activity near the sun. The range of solar 
wind speeds and the observed relations between the wind parameters 
and V suggest physical processes in the envelope other than Coulomb 
interactions. Several physical processes which might accelerate and 
heat the wind are discussed in IIIB; all of these processes occur in the 
envelope. Finally, the effects of the magnetic field and viscosity on 
the azimuthal wind speed, which is expected to be a maximum in the 
envelope, are considered in IIIC. 
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II. Structure 

A. Radial Variations 

The conditions at 1 AU and thus presumably also in the envelope are 

constantly changing, so a solar wind model can refer to any one of an 

infinite number of states at 1 AU. We shall restrict the discussion to 

some models which are designed to explain a "quiet" solar wind state, 

defined by Hundhausen (1968). It is important to understand that this 

is just one of an infinite number of states. Because n and T are related 

P 

to V, the quiet wind temperature and density are defined by specifying 

V alone, V=320 km/sec. This does not correspond to the average speed, 

most probable speed, etc; it is a state between the lowest speeds 

(s=250 km/sec) and the probable speed (f«400 km/sec). 

The density n(r), and magnetic intensity B(r), are not very sensitive 

to the thermodynamic properties of the solar wind and can be estimated 

using a 1-fluid model. Figure 2 shows n(r) and B(r) obtained by Whang 

(1970a); they are both monotonically decreasing functions. Whang's 

"quiet-time" density is consistent with the indirect measurements between 

3 

2Rq and 20 Rq and with the direct measurements (8/cm ) at 1 AU. His 
magnetic field is somewhat higher at the sun and 1 AU than the typical 
measured values (2G at the sun compared to the value 1G which is usually 
quoted, and 7.3 y at 1 AU compared to the observed value, 57), but his 
B(r) is adequate for the accuracy that we are concerned about. 

Temperatures present a greater problem, since thermodynamically the 
solar wind consists of 2- fluids and the transport parameters governing 
transfer of energy and momentum are not known for a magneto-plasma such 
as the solar wind. We shall use 2 competing models which make distinctly 
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different assumptions about the physical processes that determine the proton 
temperature T p and the electron temperature T g - the model of Hartle and 
Barnes (1970) and the model of Wolff et al. (1971). 

The proton temperature profiles predicted by the two models are 
shown in Figure 3. They both predict the same T at 1 AU, but their descrip- 
tion of the envelope is quite different. Hartle and Barnes predict a 
nearly isothermal region between 3R^and 25Rq, (4x10^ < T < 8x10 °K) because 
they assume an extended non-thermal proton heat source with characteristic 

size b=26R Q . On the other hand, Wolff et al. predict a power law decrease 

6 5 

between 3R Q and wSORq, from T = 1.2x10 °K to 2x10 °K ; they obtain the quiet 
time temperature at 1 AU by postulating an isothermal region with T e =T p = 

1.2x10 °K out to 3 R q and by i ntroc jueing viscous heating which becomes 

important at >20 Rq. The temperature rise due to viscous heating is based on 

the use of the Braginskii (1965) formula for viscosity, T)«9xl0 11 in cgs 

units. It has been argued (Parker, 1965; Hartle and Barnes, 1970) that this 

is not applicable in the region where viscous heating is said to be 

important and that the actual viscosity is negligible. 

The electron temperature at ^R^ is «1.5xl0 K and at 1 AU is «1.5xl0 K. 

Models which use the Chapman thermal conductivity predict unacceptably 

high T and conduction flux at 1 AU («3xl0 °K) (Hundhausen, 1969), but 

models which reduce the conductivity can predict a lower 

T g at 1 AU (Wolff et al. 1971) . Figure 4 show's T g (r) 

predicted by Wolff et al. (1971) on the assumption that K^K^^r/R^)^ 

9 5/2 

where K ^ is the Chapman conductivity K ec =5.5xl0 T^ and p is a free 
parameter which is taken to be p=-.728. A straight line is drawn between 
1.5xl0 6 K at 2 R q and 1.5xl0 5 at 1 AU in Figure 4 as an alternate T g (r) 


profile. 
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The relative importance of the two kinds of internal 

O 

energy - the magnetic energy density B /8rr and the thermal energy 

density nkCT^+T^) - is shown in Figure 5 which gives a plot of 
2 

Prj, = nk(T+T e )/(B /8 tt) as a function of r. This is based on the density 

profile in Figure 2, the empirical T (r) in Figure 4, the T (r) obtained 

e p 

by Hartle and Barnes (1970) and shown in Figure 3, and values of B(r) 
somewhat reduced from those in Figure 2. p increases from a very low value 
at 2 Rq (p«.l) to >1 at «50R Q . If the Wolff et al. (1971) model for T (r) 

were used (^(r) would be somewhat higher at small r and lower at higher r. 

2 

The proton f3(r), p^=nkT/ (B /8tt), shows the same general features as p^, 
but is somewhat smaller (Figure 5). 

The relation between the total internal energy density and the flow 
energy density is shown as a function of r in Figure 6. Below «15 r q the 
internal energy, which is mostly magnetic, is dominant. Above sslSR^ the 
flow energy is dominant. 

Figure 7 shows the 3 characteristic speeds of the solar wind - the 

Alfven speed, V^=B//4'np , the sound speed, Vg=V < 7p7p, and the perpendicular 

? 2 1/2 

magnetoacoustic speed, V^=(V^+Vg) as a function of r/R^, together with 
the bulk speed V(r) given by Hartle and Barnes (1970). The details of this 
figure depend on the particular models, but illustrates the following general 
results: 1) The flow becomes supersonic at s^R^ 2) V becomes comparable 

to V A and V M in the region between lOR^ and 20R Q ; 3) When V A ~ V M >>V S >>V ; 

4) When r>20R Q , ^V^V^Vg ; 5) Near 15R Q , 
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B. Azimuthal Variations 

The models discussed above are steady state, symmetric models. 
Generally, neither of these assumptions is valid for the solar 
wind. The high speed streams which recurred at 27 day intervals in 1962 
(Neugebauer and Snyder, 1966) clearly indicate azimuthal variations in 
the wind structure. Evidence for azimuthal gradients at 1 ATJ has also 
been found using widely separated spacecraft at «1 AU. Presumably, similar 
variations in the bulk speed also occur in the envelope. 

Models which would make it possible to extrapolate from observations 
of the azimuthal structure at 1 AU to the envelope do not exist. However, 
there are some idealized models which give insight into the development 
of the azimuthal structure. 

Carovillano and Siscoe (1969) considered a linear, polytropic model 
with sinusoidal velocity perturbations of V^_ and no variations in n and 

at an inner boundary at 201^. Since this model specifies the conditions 
near the edge of the envelope, it cannot be used to predict conditions 
there, but the general results are still significant. They predict 
sinusoidal perturbations in the density and azimuthal speed, growing 
linearly with distance from, the sun. The model predicts that and V^_ 
should be anticorrelated and that p should lead V by tt / (8Q) , where Q 
is the solar rotation rate and 4 streams are assumed. This qualitatively 
shows the observed snow-plow effect on density and the observed stream 
deflection, but the actual phase relations differ in magnitude from those 
predicted and the observed asymmetry in the V (t) profile is- not predicted. 
Nevertheless, the model indicates that the effects of stream interactions 
should be nearly linear near the edge of the envelope. In particular it 

will be easier to measure the V due to co-rotation if one is near the 

9 
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envelope, since the contribution to V from the stream interaction will 

cp 

be small there. Density perturbations will also be smaller near the 
envelope, so it will be easier to study the effect of solar boundary 
conditions on density with data from a spacecraft close to the sun, 

A generalization of the Carovillano-Siscoe model which includes 
non-linear effects was presented by Goldstein (1971). The results are 
illustrated in Figure 8. Near 1 AU, the model predicts the observed 
asymmetry on V , the strong density peak at the leading edge of the 
stream, and the observed form (Siscoe, 1971) of V . Clearly, non-linear 
effects are important at 1 AU. This point was also made by Burlaga et al. 
(1971). Closer to the sun, however, the results closely resemble those 
of the linear model. 

C. Time Variations 

It has long been known that the solar wind is highly variable in 
time, since ^ 27 day recurring patterns are usually not seen in the plasma 
parameters (e.g., see Figure 2 in Burlaga and Ogilvie, 1970a). Gosling 
and Bame (1971), analyzing Vela data from 1964 through 1967, found that 
few streams endure for more than 2 solar rotations and that the characteristic 
time for speed changes is about 3 days. Similar variations must occur in 
the envelope. 

The effect of a symmetric time variation in the temperature at 

tr r 

. 1 AU (a linear increase from 6x10 °K to 1.6x10 °K in «50 hrs. followed 
by a corresponding decrease) was studied by Burlaga et al. (1971) using the 
spherically-symmetric 1-fluid model of Hundhausen and Gentry (1969). 

They found an asymmetric response at 1 AU due to non-linear interactions 
(Figure 9). The general features of this profile are consistent with 
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observations, suggesting that non-linear, adiabatic compression may- 
be a basic characteristic of stream interactions. Similar non-linear 
effects are produced by a time varying spherically symmetric wind and by 
an inhomogeneous, stationary wind, so it is difficult to distinguish these 
two states; measurements of are very important in this regard. 

The effect of a linear increase in the solar wind speed at «. 1A 
over «50 hr was studied by Formisano and Chao (1971) using a spherically- 
symmetric model. They too found a non-linear interaction near 1 AU, 
and noted that the fast stream generates a pressure pulse which ultimately 
forms a shock pair at > 1 AU. This mechanism is not likely to produce 
shocks in the envelope. 

Non-linear interactions have important effects on the meso structure 
of the solar wind, but they are probably perturbations on the macroscale 
properties (Burlaga and Ogilvie, 1970 a,b). For example, Figure 10 
from Burlaga and Ogilvie, 1970a shows the distribution of positive and 
negative AV, the difference between consecutive 3-hour average values of 
bulk speed, for the interplanetary observations of Explorer 34; there is 
a difference between positive and negative gradients, presumably largely 
due to non-linear interactions, but the difference is relatively small. 
Similarly, the average wind speed, is probably not greatly changed by 
non-linear interactions, and the effect on the distributions of n and T 
is probably a skewing (Burlaga and Ogilvie, 1970 a,b) which does not 
greatly affect the long term averages. It will be interesting to compare 
the shapes of the distributions of n and T at .3 AU and at 1 AU. 
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III. Physical Processes 

A. Cosmic Ray Diffusion in the Envelope . 

Relativistic protons were emitted in an explosive solar flare on 
Feb. 23, 1956, and were observed at the earth by Meyer et al. (1956). 

Lust and Simpson (1957) noted three interesting effects in these 
observations: 1) The actual travel time of the particles was only a few 

minutes longer than the direct propagation time, 2) The particles arrived 
anisotropically from the sunward direction, 3) The more energetic particles 
arrived earlier than the less energetic particles, even though they were 
all propagating at essentially the same speed c. The first effect implies 
that particles were delayed in some region between the earth and the sun. 

If the delay were due to diffusion near the earth, then the particles 
would not have arrived anisotropically, so it is argued that they were 
delayed in a region near the sun. Now, there are 2 possibilities: Either 

the particles were trapped, for example by steady magnetic fields such as 
the magnetic field bottles identified by Schatten (1970) or by discontinuity 
surfaces (Fisk and Schatten, 1971), or they were delayed by diffusing in 
a fluctuating magnetic field. Lust and Simpson (1957) showed that the 
observed dispersion of onset times (effect 3) can be explained by diffusion 
in the region between the sun and <60 Rq. They called the region the solar 
envelope. 

Although the Feb. 23, 1956, event suggested the existence of an 
envelope filled with MHD waves, it did not provide a precise estimate 
of the size of the envelope, because all but the early minutes of the event 
were dominated by the effects of interplanetary scattering, as is the case 


for most events. 
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The May 4, 1960 event (McCracken, 1962) was ideal for studying the 
envelope because there was a strong anisotropy at 1 AU throughout the 
event, indicating that there was negligible scattering in the interplanetary 
medium at that time. Thus, the interplanetary "veil" was momentarily 
raised and we could "see" the envelope directly. Burlaga (1970) presented 
a theory for propagation under such circumstances which accurately describes 

the intensity- time profile (Figure 11) . It gives for the characteristic 

21 2 

size of the envelopeand D^IO cm /sec for the diffusion coefficient. 

Duggal et al. (1971) applied the same model to the event of Nov. 18, 1968, 
and found R^=35+15 Rq and D^«3xl0^'''cm^/sec . 

These cosmic ray observations strongly suggest the existence of 
hydromagnetic waves in the envelope, since such waves are the primary 
scattering agent at 1 AU (Jokipii, 1971) . The gyroradius of a 1 Bev 
particle near 201!^ is wSxlO^km and the Alfven speed there is «150 km/sec. 
Assuming that the particles scatter most effectively with wavelengths on 

the order of the gyroradius, \«R Q , one finds for the characteristic frequency, 
-2 -1 

(U«5xl0 sec ; this is close to the peak frequency in the observed photo- 
spheric and chromospheric acoustic spectrum (Tannenbaum et al. 1969). 

Burlaga (1970) suggested these might be the waves which Barnes (1968,1969) 
has considered to be active in heating and accelerating the solar wind. 

The technique for remote sensing of the solar envelope which was just 
described can seldom be applied at 1 AU because of interplanetary 
fluctuations, but it might be very effective closer to the sun in the region 
to be explored by the Helios and MVM spacecraft (Figure 1). 
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B. Acceleration and Heating 

1. Acceleration by heating 

a). Some basic facts. Observations at 1 AU provide 2 basic 

results which are essential clues to the acceleration and 

heating mechanisms: 1) The solar wind speed ranges from 

»250 km/sec to «850 km/sec (e.g. see the distribution in 

Figure 12), and 2) There is a quantitative relation between 

T and V which seems to be valid for all parts of the solar 
P 

cycle (Burlaga and Ogilvie, 1970a; see Figure 13). The latter 
result suggest a close relation between the heating and 
accelerating processes. 

Figure 14 compares the various kinds of energy flux and 
shows how they vary with bulk speed. This is based on the 
Tp-V relation, on averages of T g and B, (since these quantities 
do not change greatly with V - Burlaga and Ogilvie, 1970a, b; 
Ness et al. 1971; Ogilvie and Scudder, 1971), on the (weak) 
n-V relation (Burlaga and Ogilvie, 1970 b), and on observations 
of the heat flux (Montgomery et al., 1968; Ogilvie et al., 1971; 
Hundhausen and Montgomery, 1971). Clearly the most significant 
energy flux is that of the flow speed, so it is essential to 
explain the high speeds that are observed. 

b. General Discussion . The problem of accelerating the 
solar wind to high speeds (up to &850 km/sec) is a very old 
one. Parker (1963) showed that high speeds would result if 
T were constant out to some distance b <50 Rq, (Figure IS). 

Scarf and Noble (1965) pointed out that high speeds are 
difficult to explain in any other way if the coronal temperature 
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£ 

at 28^ is <2x10 K, as is generally believed. Burlaga and 

Ogilvie (1970a) showed that Parker's model produces high 

temperatures at 1 AU as well as high speeds (Table I) ; the 

results are consistent with the T -V relation within a factor 

P 

of 2, up to V«400 km/ sec. 

If one grants that high speeds are due to an extended 

region (<50E^) near the sun where remains near the coronal 

value, one is faced with the problem of explaining why T p should 

remain high. Two mechanisms have been suggested: 1) Protons 

might receive energy from electrons j which carry energy very 

efficiently from the sun; 2) Protons might be heated by some 

external source such as hydromagnetic waves. 

Consider the first mechanism. If there is a very strong 

coupling between the protons and electrons out to ^501^, then 

6 

T^T e in the envelope. Since T e «10 K out to ^301^ (Figure 3), 
one could thereby get the high T p region suggested by Parker 
(1963). But what is the nature of the coupling? The Hartle- 
Sturrock (1968) model shows that such strong coupling is not 
produced by Coulomb collisions. Hundhausen (1970) argues 
that there must be some plasma process operating in the envelope 
which maintains T^s^. The nature of this process is not 
known, but Forslund (1970) has discussed some interesting 
possibilities. 

The second mechanism, extended proton heating by an 
external source, has been reconsidered by Hartle and Barnes 
(1970) in an effort to explain the T p -V relation. They show 
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that if there is an external proton heat source, 

P (r) - D (— )exp - 
p o n ' r b 

o 

then one can explain the Tp=V relation by choosing suitable values of a, 
b and D^. They also showed the following general results: a) If energy 

is deposited below the critical point, the principle result is an increase 
in V at 1 AU ; b) If the energy is deposited above the critical point, the 
principal result is an increase in 1^; c) If the energy is deposited near 
the critical point, then both and V are increased. To explain V«350 km/sec 
and the corresponding Tp«5x]C^K (Figure 13), Hartle and Barnes (1970) 
require heating out to «£5 Kq (b=26). Burlaga and Ogilvie (1970a) obtained 
a similar result from Parker's extended heating model (Table II). 

Hundhausen (1970) has criticised the introduction of an external heat 
source on the grounds that the HS model already predicts more energy flux 
than is observed and by adding another source of energy one just make 
matters worse. However, the high flux in HS is due to the high conduction 
flux; one might be able to reduce the heat flux by using a lower K g , in 
which case the total flux predicted by Hartle and Barnes would not be 
unreasonable. Wolfe et al. (1971) have criticised the introduction 

of an external heat source on the grounds that it is unnecessary since one 
can explain the observed quiet time conditions using a viscous model with 
no external heating. However, this model does not explain the higher V, 

Tp that are more commonly observed, and the importance of viscous heating 
is still open to question . 

Both mechanisms - strong coupling ( which gives T^=T e ) and an external 
heat source (which also gives T^=«T e , within a factor of 2 - Figures 3,4) 
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imply that the solar wind behaves in some respects like a single fluid 

near the sun. Thus, in either case, a single-fluid model should give a 

reasonable zeroth order description of the solar wind near the sun. 

However, the observations that T 4 T at 1 AU (Figure 13) and that the 
7 e p 

temperature is anisotropic at 1 AU (Hundhausen, 1968) imply that the 

solar wind cannot be described as a single- fluid farther from the sun. 

Indeed, Sturrock and Hartle (1966) actually predicted that ^ before 

electron measurements were published. Thus, models such as that of 

Whang (1971a, b) with one-fluid out to b«50RQ and 2-fluids beyond appear 

to be more appropriate than 1- fluid models. 

b). Heating by Wave Damping. The Hartle and Barnes (1970) model, 

which postulates an external heat source, does not specify the physical 

nature of the source. Barnes et al. (1971) suggested that the heating 

mechanism is the collisionless damping of fast mode hydromagnetic waves 

propagating from the sun into the region of increasing [3 (Figure 5) . 

Introducing a variable iostropic flux F q of such waves with frequency 

UJ o = .02 sec ^ into the 2-fluid model of Hartle and Sturrock (1968) they 

found that the T-V relation is actually "predicted" by the model (Figure 16) 

Higher speeds and temperatures are obtained by increasing a single 

parameter, F q . Most of the wave energy is deposited within 20R G (Figure 17) 

Is there enough wave energy available to give the observed speeds? 

3 2 

The model requires F o =6.5xl0 erg/cm sec at 2 Rq to produce V=390 km/sec. 

27 

This implies an efflux of «10 erg/ sec which is less than the power 

27 

required to heat the inner corona (5x10 ergs/sec) and the chromosphere 
29 

(5x10 erg/sec). 

The model predicts that B should be nearly independent of V, in 



17 - 


agreement with observations (Burlaga and Ogilvie, 1970b) and that T g 
should decrease somewhat with increasing V. However, the densities and 

thus the total energy flux are somewhat high, especially at high speeds, 
and B is too low at all speeds. Perhaps time variations and departures 
from spherical symmetry must be considered to explain these discrepancies. 
The model also predicts a high because it uses the Chapman conductivity. 

Finally, the model does not include Maxwell stresses due to the 
background magnetic field. Since the energy in the magnetic fluctuations 
is only a fraction of that in the background field, it is not clear that 
the neglect of the latter is justified particularly where B>lG at the sun 
(Whang, private communication). 

2. Acceleration by Waves in the Envelope. 

We have seen how fast waves might accelerate the solar wind by 
damping. Alazraki and Couturier (1971) and Belcher (1971) have shown 
that Alfven waves propagating from the sun might accelerate the solar wind 
directly if they are not damped. 

The Poynting flux of Alfven waves is 
2 B 2 

F = 4rrr 7 - e(V+V A ) 
w 4tt v A' 

2 2 

where e = ( 6 B) / 2B ) and 5B is the wave amplitude. This can be very 

large near the sun because is large there (Figure 7). For example, 

6 2 2 29 

at r =10 km F :«eB V.r = e 5x10 erg/ sec which is comparable to the 
o w o A o e 

27 

total flow energy at 1 AU (10 erg/sec, Figure 14) even if e is as small 
as 1/50. 

Belcher used a 1- fluid, polytrope model to study the effect of 
such wave fluxes. Two solutions, for e=0 and for e=.01, are shown in 
Figure 18. Very low speeds at 1 AU (V«150 km/ sec) are obtained with e=0. 
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but much higher speeds (V«360 km/sec) are obtained with even small 
amplitude waves (e=.01). Most of the acceleration occurs close to the 
sun (Figure 18.). 

This model neglects damping, uses the WKB approximation, and assumes 
strong coupling between electrons and protons. A more realistic 2- fluid 
model with Alfven wave damping is currently being constructed by Barnes 
and others . 

3. Acceleration by B. 

Whang (1971a) has shown that the transformation of magnetic energy 
into kinetic energy during the expansion process might cause a small 
(~17%) increase in the solar wind speed at 1 AU at quiet times if B(r G )=2G. 
Thus, including B in a model such as that of Hartle and Sturrock (1968) 
could increase V from 250 km/sec to ^300 km/sec. 

In Whang's model the acceleration is caused by the Lorentz force 


- J x B = 

c — ~ 


4rr V 


. 2 dV 
sm cp ~r~ 
■ dr 


which is everywhere positive, directed away from the sun; here cp is the 
angle between IS and the radial direction. The acceleration results from 
the fact that near the sun the smallest cp given by Whang's 661ution, 
is 3.5 . 


The process can also be viewed as one of energy conversion. There 
is a non-zero component of the Poynting flux 

2 2 

p _ VB sin cp ^ 1 

r 4tt r^V(r) * 

Since the total energy flux is constant, the loss of Poynting flux 
ultimately appears as kinetic energy of the wind. Whang finds that 
90% of the Poynting energy is deposited with 101^ (See Figure 19). 
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C. Angular Momentum 

The flux of angular momentum is related to the Maxwell stress 
by the equation 


<?V )v - e VL 

where V is the azimuthal wind speed and L is a constant. With the 
frozen field condition and the equations VxE = V-B = 0, this gives 


Thus , 


p 

r(VB - V , B ) = -r QB = const 
cp 0 r r 



where M. = V/V. . To avoid a singularity at M =1 (the Alfven critical 
A A A 

2 

point), one must set L = Qr A , where r A is the point at which 

M =1 («25ILj. Weber and Davis (1967) find that V increases to a 
A © cp 

maximum V =3.8 km/sec at 11.5R_ and then decreases to «1 km/sec 
cpmax © 

1 AU. Thus, the effect of the magnetic torque is largest in the 
envelope, at «10R G , and is very small at 1 AD. Similar results were 
found by Brandt et al. (1969). 

Weber and Davis (1970) have also considered the effect of 
mechanical stresses on V . The basic equation in this case is 


(pV r) V = r(a° -p + -^r^ 3 ) + constant, 
cp rep r 9 I + TT 


where a is the viscous stress and p is the stress due to 
rep rep 

anisotropy. Using the Chapman formula for viscosity and assuming a 

rather high proton temperature at earth, T«2xl0^K (which greatly 

5/2 

increases the effect of viscosity, since ), Weber and Davis 

(1970) found that the viscous term is dominant outside the envelope 


and can give as large as 6 km/sec at 1 AU. On the other hand, 

using a similar model which neglects p and assumes a lower T^ 
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4 

at 1 AU (Tp!«4xl0 K) , Wolff et al. (1971) found a much lower speed 
at 1 AU (V =1.8 km/sec), which is comparable to that caused by the 
Maxwell stress alone (Weber and Davis, 1967; Brandt et al. 1969). 
Their curve for V (r), is shown in Figure 20. Of course, if Parker's 
comments on viscosity are correct, the effects of viscosity are 
negligible, and the magnitude of in Figure 20 should be reduced 
by a factor of 2. 

The observations of V are still in a rather confused state. 

9 

Experimental data for at 1 AU range from -1.5 km/sec to 

10 km/sec. (Brandt and Heise, 1970; Coon, 1969; Egidi et al., 1969; 

Strong et al., 1968). Apart from the instrumental problems of 

measuring V , there are difficulties due to the fluctuations in V 

9 9 

caused by stream interactions, (Lazarus and Goldstein, 1971; 

Siscoe, 1971; Carovillano and Siscoe, 1969) hydromagnetic waves 

(Belcher and Davis 1971) and other factors not yet understood. 

Thus, the angular momentum flux can best be studied by going close 

to the envelope or into the envelope where V is largest and V is 

cp r 

smaller. 
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IV. Discussion 

Although it is necessary to use some controversial models to 
extrapolate to conditions near the sun, it is clear that the region within 
25 Rq (the solar envelope) is very different from the interplanetary 
medium and that characteristics of the solar wind at 1 AU are largely 
determined by processes in the envelope. 

The magnetic field rules in the envelope: The Alfven speed is 

larger than the sound speed; the magnetic energy density is larger than 
the thermal energy density and the flow energy density; hydromagnetic 
waves may be propagating away from the sun, possibly scattering cosmic 
rays; hydromagnetic waves will damp, thereby heating and accelerating 
the solar wind; the Lorentz force might accelerate the solar wind to some 
extent; and the magnetic field causes the plasma to corotate. Dynamical 
effects due to stream interactions are probably small near the envelope. 

The experimental study of the envelope is of considerable importance, 
but it will remain inaccessible to direct probes for many years. In the 
meantime, the Helios and MVM spacecraft will approach the edge of the 
envelope and it will be possible to study of the envelope remotely by means 
of cosmic rays, MHD, waves etc., free of complications due to interplanetary 


effects. 



TABLE I 


b ( x )0 6 km ) 

V(km /sec ) 

T ( x I0 4 °K ) 

T (OBSERVED )(xl0 4 °K) 

5.4 

260 

0.6 

1.2 

8 

320 

1.2 

3.6 

20 

410 

5.0 

8.5 


40 


460 


14.0 


3.0 
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FIGURE CAPTIONS 


Figure 1 

Figure 2 
Figure 3 


Figure 4 
Figure 5 
Figure 6 
Figure 7 

Figure 8 

Figure 9 

Figure 10 

Figure 11 

Figure 12 

Figure 13 


Solar envelope. This shows the size and location of the 
envelope in relation to other parts of the sun's "atmosphere". 
"Quiet- time" model of n(r) and B(r). 

"Quiet-time" proton temperature profiles. Two models are 
shown. One (HB) is a non-viscous model with proton heating 
out to Rq. The other (WB) is a viscous model with T=T e 
out to 3R q . 

Electron temperature profiles. 

P(r) 

Flow energy density and internal energy density versus r. 
Characteristic speeds in the solar wind. The numbers are 
estimates which may be in error by a factor of 2. 

Stream interaction in a steady wind due to a sinusoidal 
variation in V r (cp) at IORq. 

Stream interaction due to the time varying, spherically 
symmetric heat source at .1 AU shown in the insert. 
Distribution of positive and negative AV the difference 
between consecutive 3-hour averages of the bulk speed. 

May 4, 1960, cosmic ray event. Observations are shown as 
points. The solid line is a theoretical curve. 

Distribution of speeds. This particular distribution of 
speeds is based on Explorer 34 data for the period June - 
December, 1967. 

T-V relation. Explorer 34 data for 1967 (open circles) show 
P 

a relation between the proton temperature arid speed. Data 
from other parts of the solar cycle fall on the same line, 



Figure 14 
Figure 15 


Figure 16 


Figure 17 
Figure 18 


Figure 19 


Figure 20 


indicating that the T p "V relation is a very general property 
of the solar wind. The 2-fluid model of Hartle and Barnes 
gives a point which is close to the T^-V curve. Electron 
observations suggest that is nearly independent of V, as 
indicated. 

Relative sizes of various energy fluxes versus V. 

The observed range of wind speeds can be explained by 
extended heating region. The model shown here assumes an 
isothermal region of, radius b followed by a region where 
the wind cools adiabatically. 

The lp-V relation predicted by damping of a variable flux 
F q of fast-mode waves, compared with the observed T-V 
relation. 

Most fast-mode wave energy is deposited in the envelope. 

V(r) predicted by a polytrope model with (e=0) Alfven waves. 
Alfven waves can cause significant acceleration in the 
envelope. 

Most of the Poynting energy due to a steady 1$ is deposited 
near the sun. 

V (r) , assuming both magnetic stresses and viscosity. 
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The amplification coefficient A and the solar wind velocity V as functions of r for a 

fixed value of T and two values of e (0. and .01). 
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